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ABSTRACT 

We present timing data spanning 6.4 yr for the young and energetic PSR J0205+6449, in the 
supernova remnant 3C 58. Data were obtained with the Rossi X-ray Timing Explorer, the Jodrell 
Bank Observatory and the Green Bank Telescope. We present phase-coherent timing analyses 
showing timing noise and two spin-up glitches with fractional frequency increases of ~ 3.4 x 1CP 7 
near MJD 52555, and ~ 3.8 x 1CT 6 between MJDs 52777 and 53062. These glitches are unusually 
large if the pulsar was created in the historical supernova in 1181 as has been suggested. For 
the X-ray timing we developed a new unbinned maximum-likelihood method for determining 
pulse arrival times which performs significantly better than the traditional binned techniques. In 
addition, we present an X-ray pulse profile analysis of four years of RXTE data showing that 
the pulsar is detected up to ~40keV. We also present the first measurement of the phase offset 
between the radio and X-ray pulse for this source, showing that the radio pulse leads the X-ray 
pulse by 4> = 0.10 ± 0.01 in phase. We compile all known measurements of the phase offsets 
between radio and X-ray and radio and 7-ray pulses for X-ray and 7-ray pulsars. We show that 
there is no relationship between pulse period and phase offset, supported by our measurement of 
the phase offset for PSR J0205+6449. 

Subject headings: pulsars: general — pulsars: individual (PSR J0205+6449) — X-rays: stars 



1. INTRODUCTION 

PSR J0205+6449 is a 65-ms rotation-powered 
pulsar residing in the center of the supernova re- 
mant 3C 58. The pulsar was discovered in a 
2002 Chandra X-r ay Observatory (CX O) observa- 
tion of the source (Murray et al.ll2002) and subse- 



quen tly detected as a radio pulsar (|Camilo et al 
2002). It is one of the most energetic pulsars 



in the Galaxy, with a spin-down luminosity of 
fi~ 2.7 x 10 37 ergs- 1 . 

The possible association between the 3C 58 pul- 
sar/pulsar wind nebula complex and the histori- 
cal supernova SN 1181 has been a matter of de- 
bate. 3C 58 is the only known source position- 
ally coincident and energetically compatible with 
the historic al supernova, strongly sug gesting an 
association ( Stephenson fe Greenll2002[ ) . However, 
in recent years, evidence has been mounting that 
the association may be spurious and that the true 
age of the pulsar may be closer to its character- 
istic age, r c = v/2v = 5.4 kyr (where v is the 
pulsar's rotation frequency and v its derivative), 
rather than the implied historic age of 828 yr. If 
the source is 828yr-old, the size of the remnant 
implies a large expansion velocity that is hard to 
reconcile wit h the measu red velocities of the opti- 
cal filaments ( Fesen|[l983T) a nd the expansion speed 



of the synchrotro n bubble (IBietenholz et al.ll2001 



Bietenholz|[2006h . The characteristic age estimate 
assumes that the pulsar was born spinning rapidly 
(i>q 3> v) and that the temporal spin evolution 
of the pulsar has proceeded according to a sim- 
ple magnetic dipole spin-down model. These as- 
sumptions are known to fail in some cases, such 
as for PSR J1811 — 1925 in the supernova remnant 
Gil. 2— 0.3, for which the pulsar's characteristic 
age appears to be a f actor of ~15 grea t er than that 



of th e remnant (e.g. Torii et al. 1999; iKaspi et al 



2001 



Tarn fc Roberts! 120031 ). The two ages can 
be reconciled if PSR J1811 — 1925 was born spin- 
ning slowly, with v$ ~ v. Similarly, the pulsar 
J0538+2817 has a characteristic age of 620 kyr, 
but a well established kinematic age of ~40kyr, 
i mplying a long initial spin period of ~ 138 ms 
(jKramer et al.ll2003t INk et al.ll2007h . Likewise, if 
PSR J0205+6449 was born spinning with a pe- 
riod of ~60ms, the estimated age of the pulsar 
could be reconciled with the historical supernova 



age ( Murray et al. I [20021 ). An alternate explana- 
tion of the age disparity for PSR J0205+6449 is if 
the pulsar were born with a short spin period but 
evolved more rapidly than is typically assumed for 
magnetic dipole braking. 

In principle, it may be possible to to gain infor- 
mation about a pulsar's true age via its glitching 
behavior. Pulsars in general, and young pulsars 
in particular, exhibit two types of rotational irreg- 
ularities superposed on their secular spin down. 
Glitches are characterized by a sudden increase 
in v 1 sometimes followed by an exponential decay, 
and are often accompanied by an increase in the 
magnitude of v. Glitches are common in young 
pulsars, particularly in pulsars aged ~10-100kyr, 
and provide valuable information about the su- 
pcrfluid interiors of pulsars. The fractional sizes 
of detected glitches span 5 o rders of magnitude, 
with Av/v ~ IP" 10 - 1 0~ 5 (jjanssen fe Stappersl 



20061: iHobbs et all 120021) . The nature of glitches 
appears to change with age. Glitches in the 
youngest Crab-like pulsars (1-2 kyr) tend to be 
small in magnitude, Av/v ~ 10~ 9 — 10~ 8 , and if 
the frequency recovers, it tends to recover almost 
completely. By contrast, glitches in older Vela- 
like pulsars are typically larger, with magnitudes 
of Av/v ~ 10~ 7 — 10 -6 , and typically recover 
only a small perc entage of the frequency change 
(|Lvne et al.ll2000h . 

Young pulsars are also prone to large stochas- 
tic variations in their spin-down rates, known as 
timing noise. These noise processes are observed 
as long-term trends in timing residuals after the 
removal of deterministic spin-down effects. Tim- 
ing noise tends to be the most severe in young 
pulsar s and a correlation with v is well estab- 
lished dCordes fc Downs] 1985 ; Arzoumanian et al 
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1994 ; Urama et alj 2006f ) . The low- frequency and 
broad-band nature of timing noise leads to dif- 
ficulties in its analysis. Moreover, the physical 
causes of timing noise are poorly understood. It 
could be that fluctuations in the superfluid interi- 
ors of neutron stars cause a torque on the pulsar 
(|Cordes fc Greensteinlll98ll ). or that interactions 
between t he pulsar and the magnetosphere impart 
a torque ( Chens! 1987 ). to name two possibilities. 

In this paper, we present three coherent tim- 
ing solutions for PSR J0205+6449 spanning 6.4 yr 
with data obtained from the Green Bank Tele- 
scope (GBT), the Jodrell Bank Observatory 
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(JBO) and the Rossi X-ray Timing Explorer 
(RXTE). We discuss two large glitches and tim- 
ing noise found in these data. We also present an 
analysis of the high energy emission of the pul- 
sar from the RXTE data by examining the X-ray 
pulse profile from 2 - 40keV and the phase offset 
between the radio and X-ray pulses. 

2. OBSERVATIONS AND ANALYSIS 

X-ray observations were taken with RXTE; ra- 
dio observations were taken with the GBT and 
JBO, as detailed below. The data were unevenly 
spaced throughout 6.9 yr, as shown in Figure [TJ 
and include a 202-day gap (between observations 
taken in RXTE Cycle 6 and Cycle 7) and a 
287-day gap (corresponding to Cycle 8). GBT 
data are shown with short lines, RXTE data are 
shown with medium-length lines, and JBO data 
are shown with long lines. The first set of obser- 
vations from RXTE Cycle 6 are only considered 
for the pulse profile analysis (Q and are not in- 
cluded in the 6.4 yr of data for our timing analysis 

(SO. 



2.1. RXTE Observations and Analysis 

Observations of PSR J0205+6449 were made 
using the Propor tional Counter Array (PCA; 
Jahoda et aDl2006l ) on board RXTE. The PCA is 
an array of five collimated Xenon/methane multi- 
anode proportional counter units (PCUs) operat- 
ing in the 2 - 60keV range, with a total effective 
area of approximately 6500 cm 2 and a field of view 
of - 1° FWHM. 

Data were obtained during RXTE observing 
Cycles 7, 9, and 10, spanning a period of 4yr from 
MJD 52327 to 53813, with a long gap (from MJD 
52752 to 53036) corresponding to RXTE observ- 
ing Cycle 8. In addition, four observations were 
taken in observing Cycle 6 (17.1 hr exposure over 
MJD 52138 - 52141). Because of the long data 
gap between the Cycle 6 data and the subsequent 
observations in Cycle 7 (202 days), they could 
not be unambiguously phase connected, render- 
ing these observations of limited interest for our 
timing analysis. They were considered for a pulse 
profile analysis discussed in Sj4] The frequency for 
each Cycle 6 observation was determined by a peri- 
odogram analysis, which was then used for folding 
the data. All of the RXTE data were collected 



in "GoodXenon" mode, which records the arrival 
time (with l-/zs resolution) and energy (256 chan- 
nel resolution) of every unrejected event. Typi- 
cally, 3 PCUs were operational during an obser- 
vation. For our timing analysis, we used only the 
first layer of each operational PCU in the energy 
range 2 - 18keV, as this maximizes the signal- 
to-noise ratio of individual observations for this 
source. 

Observations were downloaded from the HEASARC 
archival! and photon arrival times were converted 
to barycentric dynamical time (TDB) at the solar 
system barycenter using the J2000 source posi- 
tion determined using data from CXO, RA = 
02 h 05 m 37!92 ± Of 02, Dec = 64° 49' 42.8" ± 0.72" 



(jSlane et al.ll2002D and the JPL DE200 solar sys- 
tem ephemeris with the FITS tool 'faxbary'. 

We noted that two observations directly follow- 
ing the leap second occurring on 2006 January 1 
had incorrect clock corrections, confirmed by the 
RXTE team (C. Markwardt, private communica- 
tion). This was fixed by adding a Is time jump 
to each pulse time-of-arrival obtained from these 
observations. 

2.1.1. A New X-ray Timing Technique 

The X-ray pulse profile as measured in each ob- 
serving session by RXTE has two narrow compo- 
nents separated by approximately one-half of a ro- 
tation of the pulsar (See Fig. [2]). Such sharp pulse 
features are usually very helpful for timing analy- 
ses, however, the "signal-to-noise" ratio (here de- 
fined as the ratio of pulsed source counts to all 
other counts) for the RXTE observations was very 
low, with typical values near 6xl0 -3 . These lev- 
els are so low that for some observations the pulses 
were not easily visible in the binned pulse profile 
plots. Therefore, in order to determine the pulse 
times-of-arrival (TOAs) from the X-ray data, we 
used a new and simple maximum-likelihood (ML) 
technique. This method avoids the information 
loss inherent in binning event-based (i.e. photon) 
pulse profiles. The inputs to the method are an 
accurate model pulse profile /($) of intensity, in- 
cluding both signal and the average noise level or 
background, as a function of rotational phase <& 
(where < $ < 1), and the computed rotational 
phases 4>i of the N events (where 1 < i < N) 
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from the observation according to the best timing 
model of the pulsar. We normalize /($) so that 
it has unit area and can be treated as a probabil- 
ity density function for the individual event arrival 
times. 

If we assume an arrival offset t/P (where P 
is the spin period) of rotational phase (where 
< t/P < 1), we can compute a probability or 
likelihood that our data has that offset from the 
template using Prob(t/P) = n£Li !{<t>i ~ t / p )- 
If we compute a large grid of 1000-10000 evenly- 
spaced probabilities for offsets t/P between and 
1, the resulting distribution describes the proba- 
bility density for the average pulse arrival time. 
When normalized, we can determine the most 
likely arrival time, or by suitably integrating the 
distribution, the median arrival time and error es- 
timates for the arrival time. 

To test how well the ML technique determines 
TOAs compared to the traditional binned Fourier- 
base d methods ty pically used in radio pulsar tim- 
ing (|Tavlorl[l99i l. we simulated ~100,000 high- 
energy observations of pulsars with both a pulse 
profile like that of PSR J0205+6449 (i.e. two nar- 
row Gaussian peaks separated by nearly half a ro- 
tation, Fig. [2]) and a wide Gaussian profile with 
FWHM=0.2, each with a wide variety of signal- 
to-noise ratios, total photon counts, and numbers 
of bins in the profiles used to measure the pulse 
phase for the traditional Fourier method. 

The simulations showed that: 1. The tradi- 
tional Fourier-based technique underestimated the 
TOA errors by 10% to several hundred percent de- 
pending on the signal-to-noise ratio of the profiles 
as well as the number of bins used in the pulse 
profile. In contrast, the average errors estimated 
by the ML technique were almost always within 
5%, and were typically within 1—2%, of the true 
TOA errors as determined from the statistics of 
the simulations. 2. The measured TOA error 
distributions were narrower for the ML technique 
than for the Fourier technique (i.e. the TOAs were 
more accurate) by several percent typically, but 
by up to many tens of percent for low signal-to- 
noise ratio cases where the ratio of pulsed counts 
to background counts is below ~ 0.01 — 0.05. 3. 
In very low signal-to-noise cases (i.e. where the 
ratio of pulsed counts to background counts is 
~ 0.001 — 0.01), the Fourier method determined 
completely wrong TOAs (i.e. where the measured 



TOA differed from the true TOA by many times 
the estimated TOA error) much more often than 
the ML technique, typically by factors of 2-5 times. 

In summary, the ML technique determines 
slightly more accurate TOAs, with much better 
error estimates, over a wider range of signal-to- 
noise ratios, than the traditional binned Fourier 
technique of TOA determination. We recommend 
that it be used in all X-ray and 7-ray pulsar tim- 
ing applications, particularly in low signal-to- noise 
ratio situations. 

2.1.2. Application to PSR J0205+6U9 

For the ML-timing of PSR J0205+6449, we 
used a two-Gaussian profile model or template 
/($) based on a high signal-to-noise ratio pulse 
profile from four months of RXTE data from early 
2004 (Cycle 9). The timing model used for the 
profile included several frequency derivatives such 
that no timing noise was apparent in the timing 
residuals. The two-Gaussian plus DC component 
model was then fit to the data. The parameters 
of the resulting model are shown in Table [TJ The 
Gaussian template with 1000 phase bins is shown 
in Figure [2j 

For each RXTE observation, we folded the 
X-ray data using the best predicted spin period 
for that day w ith the software package PRESTO 
(jRansoml 120011 ). but we allowed the software to 
search in a narrow range for the best pulsation 
period. This technique was necessary to optimize 
the signal-to-noise ratio of the folded profiles be- 
cause of the large levels of timing noise from the 
pulsar. We then used the refined spin period to 
determine pulse phases <f>i for each of the X-rays. 

To determine the ML-derived TOAs, we used a 
grid of 1000 evenly spaced phase offsets t/P over 
a full rotation of the pulsar. While this is a fairly 
computationally expensive task when the number 
of events is large (as is the case for RXTE with 
its high background rate), the resulting likelihood 
distribution was typically well behaved (i.e. uni- 
modal and nearly symmetric) and produced ex- 
cellent TOAs and error estimates. Typical uncer- 
tainties were between 450— 750 ns for each TOA. 

The phase offset, as determined by the median 
of the resulting likelihood distribution, was multi- 
plied by the current pulse period and added to the 
fiducial epoch for each observation, in our case the 
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first X-ray recorded, to make a TOA. Approximate 
one-sigma error estimates for each TOA were de- 
termined by integrating the likelihood distribution 
in each direction until 0.8413 of the total likelihood 
was accounted for. For each RXTE observation, 
we typically determined 2—3 TOAs. The result- 
ing TOAs were fitted to a timing model (see §3.1 
and §3.2) using the pulsar timing software package 
TEMPC0. 

2.2. Green Bank Telescope Observations 

Observations with the GBT were made at ei- 
ther 820 or 1400 MHz from MJD 52327 - 52776 us- 
ing the Berkeley-Ca ltech Pulsar Machine (BCPM 



to UTC using informati on from the GPS. F urther 



Back er et all 11997 ). The BCPM is an ana- 



log/digital filterbank which samples each of 2x96 
channels using 4 bits at flexible sampling rates 
and channel bandwidths. We recorded data using 
134 MHz of bandwidth and 50 /is samples for the 
1400 MHz observations or 48 MHz of bandwidth 
and 72 /is samples for the 820 MHz observations. 
Typical integrations times lasted between 3—5 hrs. 
We de-dispersed at the know n dispersion measure 



(DM) of 140.7±0.3pccm" 3 (jCamilo et al. 2002) 



and f olded all of the data using PRESTO (| Ransom 
l200l[ ). and then extracted 2—3 TOAs per obser- 
vation by correlating in the frequency domain the 
folded 64-bin pulse profile with a Gaussian of frac- 
tional width 0.04 in phase. The typical precision of 
the TOAs was 200-400 /is. TOAs were corrected 
to the UTC timescale using data from Global Po- 
sitioning System (GPS) satellites. 

2.3. Jodrell Bank Observatory Observa- 
tions 

Observations were made at JBO every 3—5 days 
between MJD 53725 and 54666 using the 76-m 
Lovell telescope of the University of Manchester at 
a frequency of 1.4 GHz. Each observation typically 
lasted 3 hr, divided into 1-minute sub-integrations. 
The data were de-dispersed at the known DM in 
hardware and folded on-line. The profiles from 
these, sampled in intervals of 164.3 /is, were added 
in polarization pairs and then combined to provide 
a single total-intensity profile. This was then con- 
volved with a template derived from a single high 
signal-to-noise-ratio 400-bin profile at the same 
frequency to yield a TOA. TOAs were corrected 



details can be found in iHobbs et alj (|2004f ). The 
typical precision of the TOAs was 200— 700 /is. 

3. TIMING ANALYSIS 

3.1. Phase-Coherent Timing Analysis 

The most accurate method used to determine 
spin parameters is phase-coherent timing, that is 
comparing TOAs with a model ephemeris and 
accounting for each rota tion of the pulsar, as 
descr ibed elsewhere (e.g. Lvne fc Graham-Smithl 
2005). A single phase-coherent timing solution 
spanning all 6.4 yr of data proved impossible for 
these data owing to the 287-day gap in timing ob- 
servations as well as two large glitches. As such, 
we present three coherent timing solutions, sum- 
marized in Table [U We verified that changing the 
pulsar position by 3c did not significantly change 
the fitted parameters, thus, the pul sar position 
was h eld fixed at that determined by ISlane et al 
(|2002t ) from CXO data. 
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GBT and RXTE data were fitted together, re- 
sulting in a timing solution spanning six months 
(MJD 52327 - 52538), with GBT (dots) and 
RXTE (crosses) as shown in Figure [3] Timing 
residuals are shown in the Figure, with the top 
panel showing residuals with v, z>, and v fitted. 
Significant timing noise remains in the data and 
can be fitted with six frequency derivatives, shown 
in the bottom panel of Figure [31 

Phase coherence was lost after MJD 52538, as 
the result of a glitch (see §3.2j) . A second coher- 
ent timing solution using GBT and RXTE data 
spans seven months (MJD 52571 - 52776). Tim- 
ing residuals are shown in Figure^ with is, v, and 
v fitted in the top panel. As with the previous 
timing solution, significant timing noise remains 
in the residuals, which is fitted with 5 frequency 
derivatives, shown in the bottom panel of Figured] 

Figure [5] shows residuals for the third timing so- 
lution spanning -4.4 yr (MJD 53063 - 54669). In 
addition to two years of X-ray timing observations 
(spanning MJD 53063 - 53813), on MJD 52725 
(2005 December 21), radio timing observations be- 
gan using JBO. X-ray and radio observations were 
concurrent for 88 days, after which no additional 
X-ray observations were obtained. 

The top panel of Figure [5] shows residuals with 
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v, i>, and v removed; the bottom panel shows 
residuals with 12 frequency derivatives fitted. The 
timing noise in this 4.4-yr period is so large that 
it cannot be fully described by a 12 degree poly- 
nomial (the largest allowed with current machine 
precision). In addition, some of the timing noise 
seen in these data is likely attri butable to unmod- 
eled glitch recovery (|Lvnelll996h . In fact, the mea- 
sured value of v is significantly different from that 
predicted from the previous timing solution. The 
difference between the predicted and measured v is 
too large (~6x 10~ 5 Hz) to be explained by timing 
noise, indicating that a glitch probably occurred 
during the 287-day gap in the data. 

3.2. Glitches 

In general, the frequency evolution of a pulsar 
following a glitch can be characterized by 

v(t) = m(t) + Ai/[1 - Q(l - e-^ Td )} + Az> p i, (1) 

where vo(t) is extrapolated from the pre-glitch 
value, Av is the initial jump in frequency, Q is 
the recovery fraction, Td is the time decay con- 
stant, and Ai> p is the permanent change in v as a 
result of the glitch. 

In order to analyze further the two glitches in- 
ferred from our coherent timing analysis, we per- 
formed a partially coherent timing analysis over 
short time intervals (on average ~40days), fit- 
ting only for v and v and choosing the length of 
each data subset such that the phase residuals are 
Gaussian-distributed (i.e. 'white'). The results, 
with the average v from the inter-glitch period re- 
moved, are shown plotted in the top panel of Fig- 
ure [6] Two glitches as well as timing noise are 
clearly present. We also show measurements of v 
for the same intervals in the bottom panel of Fig- 
ure [B] Again, two glitches are apparent, as is the 
significant timing noise in the data. 

To measure accurately the size of each glitch 
while minimizing the contaminating effect of long- 
term timing noise, we took v measurements span- 
ning only ~200 days before and after each glitch to 
measure the fractional increase in spin frequency. 
Figures [7] and [8] show pre- and post-glitch v mea- 
surements with the post-glitch slope subtracted. 

We observed a frequency increase between 
MJDs 52538 and 52571 (Fig. E}, corresponding 
to the loss of phase coherence discussed in §3.11 



This v increase corresponds to a glitch of frac- 
tional magnitude Avjv = (3.4 ± 1.1) x 1CT 7 
(see also Ransom et al. 2004). The change in 
v over the glitch is not statistically significant 
when measured as the slope of frequency mea- 
surements, as shown in the Figure. However, the 
fractional change in v from the individual v mea- 
surements before and after the glitch is significant, 
with Au/u ~ 0.005 ± 0.001. No short-term post- 
glitch relaxation is detected; however, because 
of the sparse sampling, it cannot be precluded. 
Neither can a long-term post-glitch relaxation be 
distinguished from a simple change in v because 
a decay, if present, was interrupted by a second 
glitch. 

We observed a second frequency increase be- 
tween MJDs 52776 and 53063 as shown in Fig- 
ure [H Because no timing data were taken during 
the 287-day period where the frequency jump oc- 
curred, it is not possible to differentiate between a 
single glitch and two (or more) smaller glitches. 
However, because there is no clear evidence of 
more than one glitch, we interpret the frequency 
increase as a single glitch of fractional magnitude 
Avjv = (3.8 ± 0.4) x 10" 6 . The frequency deriva- 
tive, as measured from the slope before and af- 
ter the glitch (see Fig. [8]) also changed signifi- 
cantly, with a fractional magnitude of Av/i> = 
0.012 ±0.001. 

Another possible description of this glitch is 
a large increase in v followed by an exponen- 
tial recovery (also suggested for this glitch by 
W. Hermsen, personal communication). Fitting a 
glitch model including an exponential decay to all 
the frequency measurements before and after the 
glitch, we found models which span a wide range 
of possible glitch parameters, with a reduced x 2 of 
~ 3.6 — 3.9 for 25 degrees of freedom. We found fits 
of roughly equal probability for all possible glitch 
epochs between the two bounding coherent timing 
solutions, thus we present a range of glitch pa- 
rameters corresponding to the date limits of MJD 
52777 and MJD 53062. A typical fit and residu- 
als are shown in Figure [9l The fractional magni- 
tude of the glitch from these models ranges from 
Avjv ~ (3.6 — 7.2) x 10~ 6 , while the recovery frac- 
tion spans Q ~ 0.66-0.88 and r d ~ 280-295. The 
long-term change in frequency derivative, Av /£>, is 
always negative, i.e. in the opposite direction from 
that expected in a typical glitch, and is approxi- 
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mately equal to ^—0.0046 for all fitted models. 
This could be attributed to unusual glitch recov- 
ery, timing noise, or a combination thereof. How- 
ever, this effect is possibly an artifact of the fitting 
procedure, where it is assumed that the pre-glitch 
v is not itself recovering from the previous glitch. 

Because of the unusual behavior of the long- 
term Av after the glitch, and the possibility that 
this behavior is not a direct result of the glitch, 
we performed a third fit to the frequency data, 
this time excluding frequency measurements af- 
ter MJD 53700 (after which the unusual Av dom- 
inates), and fixing the change in v to be zero. 
Fitting an exponential glitch recovery model to 
this subset of data provides a much better fit, 
as shown for a sample glitch epoch in Figure 1101 
Again, the glitch may have occurred at any time 
between MJD 52777 and 53062. For this fit, the 
range of possible fractional increases in frequency 
is Av/v = (3.4 — 5.3) x 10~ 6 , while the recovery 
fraction is Q = 0.5 — 0.67 and the recovery time 
scale is = 379 — 383 days. 

4. X-RAY PROFILE ANALYSIS 

For each RXTE observation in observing Cycles 
7, 9, and 10, we created a phase-resolved spectrum 
with 64 phase bins across the profile, using the 
Ftool 'fasebin' and the partially coherent timing 
ephemerides described above ( §3.2|) . We also cre- 
ated phase-resolved spectra as described above for 
the four observations taken during RXTE Cycle 
6. These observations could not be unambiguously 
phase-connected with the later data because of a 
202-day gap between these and the subsequent Cy- 
cle 7 observations. Instead, we folded the observa- 
tions at a frequency determined by a periodogram 
analysis, v = 15.22466(2) Hz with a frequency 



derivative of v 



-4.49(1) x 10 



evaluated 



at the reference epoch MJD 52139. 5. Because of 
the h ard spectrum of the source ( Ransom et al.1 
20041 ). we selected photons from all three layers 
of the Xenon detectors, and from all operational 
PCUs and included all photons within the energy 
range 2 - 60keV. In addition, we repeated the 
entire process using the first Xenon layer for the 
energy range 2 - 18keV, where the majority of 
the softer source counts reside. We also created 
phase-averaged spectra which were used to build 
response matrices for each PCU. We then used 



XSPECd to create a pulse profile in count rate for 
each PCU for each observation. We then combined 
data from all PCUs for each observation, dividing 
by a factor to account for the amount of time each 
PCU was on, resulting in a pulse profile in units 
of count rate per PCU. 

We aligned the 2 - 18keV profile for each ob- 
servation by cross-correlating with the 2-Gaussian 
template (Fig. \2§ and summed profiles from each 
RXTE cycle together, as shown in Figure QT] 
There is no evidence for evolution in the profile 
over ~4yr. 

We also created profiles for different energy 
bands: 2 - lOkeV, 10 - 18keV, and 18 - 40keV 
and summed profiles from RXTE cycles 7, 9, and 
10 (where phase-coherent timing solutions exist) 
into a single profile for each energy band, shown 
in Figure [12] We used the phase offset determined 
from cross correlating the 2 - 18keV profile to 
align the profiles from all the energy ranges. This 
ensures that the correct phase offset is applied the 
18 - 40keV profiles, where the signal-to-noise ra- 
tio of individual profiles is very poor, and a noise 
spike can be mistaken by the cross correlation al- 
gorithm for the sharp peak of the profile. We note 
that using the phase offset determined for a dif- 
ferent energy range profile could be problematic 
in building added profiles if there is signifcant en- 
ergy evolution across the energy range of interest. 
This does not appear to be the case for this pulsar 
for the relevant energy range. The pulsar is visible 
up to ~ 40 keV and the pulse shape does not vary 
significantly with increasing energy. 

There is some possible unusual structure in the 
off-pulse region of the profiles, particularly for the 
energy range 18 - 40keV. We examined this pu- 
tative structure by performing several additional 
analyses. First, we produced summed profiles us- 
ing a different method than described above. We 
extracted events and created a time series for each 
energy range, for each individual observation. We 
folded each time series with the local ephemeris 
and aligned and summed the profiles several times 
using different template profiles. This analysis 
produced pulse profiles that are not significantly 
different from those shown in Figure 1121 

We further analyzed the off-pulse region of 
the profiles by applying a \ 2 test, a Z 2 test 



3 See http://xspec.gsfc.nasa.gov, ver. 11.3.1. 
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(IBuccheri et al ll983h . and an H test |De Jager et al 
1989J). The reduced x 2 values for the 10 - 18keV 
and 18 - 40keV are less than 1, while the % 2 for 
the off-pulse region of the 2 - lOkeV profile is 55 
for 48 degrees of freedom (where the probability 
of this x 2 or higer occurring by chance is 21%). 
The Z 2 test for 1, 2, 4, and 8 harmonics, and the 
H test applied to the off-pulse region of all three 
profiles resulted in the null hypothesis. There- 
fore, the off-pulse regions of the profiles are not 
statistically different from a DC offset. 

5. PHASE OFFSET BETWEEN THE 
RADIO AND X-RAY PULSES 

Precise measurement of the phase lag between 
the radio and X-ray pulses is important for under- 
standing the pulse emission mechanism. Emission 
from rotation-powered puls ars is thought to arise 
from either a polar cap (e.g. Daughertv k. Hardingl 



1982 ), or in magnetospheric o uter gaps (e.g. 



Cheng et ai1ll986t (Romanilll996l ). Absolute tim 



ing for several different energy ranges can place 
constraints on the shape of the outer gap and 
the height in the magnetosphere where ra diation 
is generated (jRomani k, Yadigaroglul fl995h . The 
radio pulse profile of PSR J0205+ 6449 is single 
peaked with a width of ^0.05P ( Camilo et al 



2002); the X-ray profile is double peaked, with a 
peak-to-peak separation of ^0.5. We made two 
independent measurements of the phase offset be- 
tween the radio pulse and the main X-ray pulse 
for PSR J0205+6449, by finding the offset be- 
tween the GBT and RXTE data, and the JBO 
and RXTE data. 

We used the first timing solution (GBT and 
RXTE data prior to glitch 1) spanning MJDs 52327 
to 52538 to make three independent measurements 
of the phase offset. We split our timing solution 
into three subsets, fitting only for v and v. Ra- 
dio TOAs were sh ifted to infinite freq uency using 



the nominal DM IjCamilo et all 120021 ). The data 



subsets were chosen such that there was good 
overlap between the sparsely sampled X-ray and 
radio data, and that each solution had Gaussian- 
distributed residuals. The weighted average value 
of these three measurements is 6.72 ± 0.66 ms, 
or 0.102 ± 0.010 in phase (radio leading). This 
analysis could not be repeated for the post-glitch 
GBT / RXTE timing solution because of the even 



more sparsely sampled data, which consist of sig- 
nificant lags between most radio and X-ray ob- 
servations, leading to poorly constrained values of 
the phase offset for short sections of the data. 

We obtained a second measurement of the 
phase offset by simultaneously fitting overlapping 
timing observations from RXTE and JBO, span- 
ning 88 days from MJDs 53725 to 53813. We first 
obtained a phase-coherent timing solution from 
the well sampled radio data, and then added the 
overlapping, more sparsely sampled X-ray TOAs. 
Again, radio TOAs were shifted to infinite fre- 
quency using the nominal DM value. We split 
the data into three subsets, fitting only for v and 
v and ensuring that each subset had Gaussian 
residuals. The JBO TOAs were created using a 
different fiducial point than the GBT and RXTE 
profiles, so each TOA was shifted by a constant to 
adjust for the difference between the two fiducial 
points used. The weighted average of these offset 
measurements is 5.55 ± 0.66 ms, corresponding to 
a phase offset of 0.085 ± 0.010. The difference 
between this measurement and that made with 
RXTE and GBT data is 1.17 ± 0.94ms, i.e., in 
agreement within 1.25a. 

We further verified the phase offset by extract- 
ing a single TOA from an archival CXO observa- 
tion from February 2002 (Observation ID 2756). 
The offset measured from the CXO TOA to the 
GBT radio TOAs agrees with our GBT / RXTE 
offset within 1.2 a. 

The uncertainty in the phase offset is domi- 
nated by the uncertainty in the DM. The DM 
was measured to be 140.7±0.3pc cm~ 3 using 800 
and 1375 MH z GBT data in 2002, as reported in 
Camilo et all (|2002l ). Thus the JBO determina- 
tion of the phase lag, made ~3.5yr after the mea- 
surement of DM, could be affected by short or 
long-term changes in DM. We therefore report the 
first measurement of the phase offset made with 
GBT and RXTE, of <j> = 0.102 ± 0.010. 

6. DISCUSSION 

6.1. Timing Noise, Glitches, and the Age 
of PSR J0205+6449 

We show evidence for two large glitches in 6.4 yr 
of timing data of PSR J0205+6449. The fractional 
magnitudes of these glitches (Av/v ~ 10~ 7 — 10~ 6 ) 
are typical of pulsars with characteristic ages of 
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5— lOkyr, such as the frequent large glitcher PSR 
J0537-6910 (r c = 4.9 kyr) or the Vela pulsar (r c = 
llkyr). The youngest pulsars such as the Crab 
pulsar (955 yr), PSR B0540-69 (r c = 1.7 kyr) and 
PSR J1119-6127 (r c = 1.6 kyr) typically have 
glitches with smaller fractional magnitudes rang- 
ing from Au/u ~ 10~ 9 — 10~ 8 . Perhaps this indi- 
cates that the pulsar age is closer to its character- 
istic age of t c ~ 5.4 kyr, rather than the historical 
supernova age of 828 yr. 

If the pulsar was born in the historical super- 
nova event 828 yr ago, these glitches are unusu- 
ally large. Glitches may be rel ated to pulsar age 
via n eutron star temperature (jMcKenna fc Lyne 
1990). If this is the case, the large glitches ob- 
served here could be related to the ver y low mea- 
sured temperature of PSR J0205+6449 (|Slane et al 



2004), rather than its chronological age. The rea- 
son for the exceptionally cool surface tempera- 
ture of this neutron star is still a mystery, though 
may be explained wit h a large neutron star mass 
(jYakovlev et al.ll2002l ). On the other hand, large 
glitches have been observed in the hot surface 
tem perature Anomalous X-ray Pulsar s (AXPs, 
e.g. iKaspi et alj l2000t iDib et alj 120071) ■ If the 
mechanism behind rotation-powered pulsar and 
magnetar glitches is similar, the neutron star sur- 
face temperature may not be the primary factor 
in determining the size of glitches. 

Young pulsars emit large amounts of energy 
as they spin down, providing for easy measure- 
ment of i>, and occasionally hi gher order fre- 
quency derivatives (v, v; e.g. iLvne et al 1993t 
Livingstone et al ] l2005h . A measurement of v pro- 
vides an estimate of the 'braking' index, n — 
vi)/i> 2 , giving insight into the physics underlying 
pulsar spin-down, as well as an improved esti- 
mate of the pulsar age. Both timing noise and 
glitches can prevent a measurement of n. Gener- 
ally, only the youngest pulsars, with t c < 2 kyr, 
have measurable braking indices. The exception 
is the 11 kyr-old Vela pulsar, where frequent large 
glitches prevent a phase-coherent measurement of 
the braking index, but measuring v in the after- 
math of glitches has allowed a meas urement of 
n = 1.4 ± 0.2 over ~25yrs of data (jLvne et al 
19961) . 



three phase-coherent timing solutions obtained for 
this source, ranging from n ~ 15 — 90. A partially 
coherent timing analysis (Fig. [6]) shows that v does 
not evolve linearly implying that a deterministic 
value of v and thus n cannot be measured from 
these data. However, excluding data in the imme- 
diate aftermath of the glitches and looking at the 
overall trend in v from the first three and last 10 
measurements of v from Figure [6] (bottom panel), 
the implied value is n ~ 4. Though this value is 
contaminated by timing noise and glitch recovery, 
it is suggestive that the true, underlying value of 
n may eventually be measurable with long-term 
timing. 

6.2. Absolute Timing and the Pulse Emis- 
sion Mechanism 

The observed phase difference between radio 
and high-energy pulses, as well as the pulse shape 
and peak-to-peak separation, in principal provide 
important information about the pulsar emission 
mechanism by constraining the pulse emission re- 
gion. Table [3] shows the offsets between the radio 
and X-ray pulses and the radio and 7-ray pulses 
for all known measurements to date. These com- 
piled data show that there is a large scatter in the 
measured offsets and no correlation between pulse 
period and phase offset. 

The measured offset for PSR J0205+6449 
of <fi = 0.10 ± 0.01 supports the lack of cor- 
relation between phase offset and pulse pe- 
riod. In particular, PSR J0205+6449 has a 
very similar pulse period (65.7 ms) to the pul- 
sar PSR J1420-6048 (68.2ms), while the phase 
offset of PSR J1420-6048 is = -0.35(6) 



(jRoberts et al.ll200lh . It is likely that the geome- 
try and viewing angle of each system will affect the 
measured phase offset and may ultimately explain 
the range of observed offsets. 

The radio-to-X-ray phase offset for PSR J0205+6449 
is consistent with the radio-to-7-ra y offset of 



0.08 ± 0.02 (|Abdo et alJ l2009bh . The X 



The initial goal of timing PSR J0205+6449 was 
to measure its braking index. However, the mea- 
sured value of n varies significantly among the 



ray and 7 -ray phase offsets are li kewise aligned for 
the Crab dPellizzoni et al. 2009a ) and Vela pulsars 
(e.g. lAbdo et al.ll2009dl ). while this i s not the case 
for th e young pulsar PSR B1509— 58 (jKuiper et al 



1999) , nor the millis econd pulsar PSR J0218+4232 
(lAbdo et al.ll2009al ). 



The outer gap model of lRomani fc Yadigaroglu 
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(|l995l ) predicts that the X-ray pulse should lag 
the radio pulse by 0.35—0.5 in phase and should 
appear as a single broad pulse. Neither pre- 
diction is supported by the X-ray profile of 
PSR J0205+6449. In addition, a thermal X-ray 
pulse arriving in phase with the radio pulse is pre- 
dicted, while no thermal pul sations have been d e- 
tected in PSR J0205+6449 (|Murrav et alJl2002h . 



The main pulse of PSR J0205+6449 is detected 
up to 40 kcV with the PCA on board RXTE, while 
the interpulse is visible up to ~18keV. It is one of 
an increasing number of young pulsars detected in 
the hard X-ray energy range and has also recently 
been detect ed up to ~ 3 GeV w ith the Fermi Space 
Telescope (|Abdo et all l2009bl ). Comparing the 
emission of the Crab pulsar and nebula to the 3C 
58 pulsar and nebula at higher energies may offer 
insight into the physical reasons behind the in- 
triguing differences between these two seemingly 
similar objects. 

7. CONCLUSIONS 

Multi-wavelength timing observations offer an 
excellent probe of both the temporal and emis- 
sion characteristics of young pulsars. We observed 
two large glitches that are not characteristic of the 
proposed young age of PSR J0205+6449. This is 
not conclusive however, and the timing data are 
consistent with an age of 828 yr if the pulsar was 
born spinning slowly with Pq ~ 60 ms, or an age 
of r c ~ 5.4 kyr if the pulsar was born spinning 
rapidly. Furthermore, the age of the pulsar is con- 
sistent with r > t c = 5.4 kyr if its true braking 
index n < 3, as is the case for all measured val- 
ues of n, though there is currently no evidence for 
n < 3 for PSR J0205+6449. Long-term timing 
may eventually allow for the measurement of n 
by using the incoherent meth od performed on the 
Vela pulsar (jLvne et al. I ll996h . 

We have presented the first measurement of the 
phase offset between the radio and X-ray pulses 
of PSR J0205+6449 to be 0.10 ± 0.01, which is 
consistent with the recently reported 7-ray phas e 
offset of <f> = 0.08 ± 0.02 lAbdo et alJ (l2009bh . 
PSR J0205+6449 is rare in that it has both the 
magnetospheric X-ray and 7-ray phase offset pre- 
cisely measured. Among the pulsars with mea- 
sured phase offsets, with periods ranging from 
1.56 ms to 1250 ms, there is no correlation between 



pulse period and phase offset, as shown in Ta- 
ble [3] Phase offset measurements should be im- 
portant for constraining the pulse emission region 
and significant progress is occurring at present in 
the measurement of radio-to-7-ray phase offsets 
with the ongoing detection of many radio pulsars 
as 7-ray pulsars with Fermi. 
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Fig. 1. — Distribution of observations of 
PSR J0205+6449 from RXTE, GBT, and JBO 
over 6.9 yr. GBT data are indicated with short 
lines, RXTE data with medium length lines, and 
JBO data with long lines. The three coherent tim- 
ing solutions, spanning a total of 6.4 yr, are indi- 
cated with horizontal lines along the bottom of the 
plot. 
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Fig. 2. — Best- fit two-Gaussian curve with 1000 
phase bins from maximum-likelihood method of 
determining pulse times of arrival for 2 - 18keV 
as described in 32.1.11 



Fig. 3. — Timing residuals for RXTE observing 
Cycle 7 (crosses) and GBT data (dots) spanning 
MJDs 52327 - 52538. The top panel shows resid- 
uals with v, v, and v fitted. The bottom panel 
shows residuals with an additional four frequency 
derivatives removed. 
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Fig. 4. — Timing residuals for RXTE observing 
Cycle 7 (crosses) and GBT data (dots) spanning 
MJDs 52571 - 52776. The top panel shows resid- 
uals with v, v, and v fitted. The bottom panel 
shows residuals with an additional three frequency 
derivatives removed. 
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Table 1 



X-ray Profile Template Parameters 
Parameter Value 

Flux 1 0.64515 
FWHM 1 0.02386 
Phase 1 0.0 
Flux 2 0.35485 
FWHM 2 0.05826 
Phase 2 0.50515 
DC flux 193.46 

Parameters for the two-Gaussian template profile, which is proportional to the X-ray arrival time 
probability density function used to determine RXTE times-of-arrival (see §2.1.1j) . The peak for the first 

Gaussian was explicitly placed at zero. The listed values for flux indicate the integrated area of each 
Gaussian such that the total pulsed signal "flux" from the pulsar equals one. The DC flux describes the 

integrated background level over the full pulse profile. 
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Fig. 5. — Timing residuals for RXTE observing 
Cycles 9 and 10 (crosses), and Jodrell Bank data 
(small dots), after the glitch occurring between 
MJD 52777 and 53062. The top panel shows resid- 
uals with only v, i>, and v fitted. The bottom 
panel shows residuals with v and 12 frequency 
derivatives fitted (the maximum allowed with cur- 
rent machine precision). 



Fig. 6. — Frequency and frequency derivative evo- 
lution of PSR J0205+6449 over 6.4 yr. The top 
panel shows frequency measurements from short, 
phase-coherent timing solutions with the over- 
all trend in the inter-glitch interval subtracted. 
The bottom panel shows measurements of the fre- 
quency derivative. Uncertainties are smaller than 
the plotted points. 
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Table 2 

Phase-coherent Timing Parameters for PSR J0205+6449 
First phase-coherent solution a 



Dates (Modified Julian Day) 
Dates 

Number of TOAs 

Epoch (Modified Julian Day) 

v (Hz) 

v (ltr 11 s- 2 ) 
v (1(T 21 s- 3 ) 

RMS residuals with i) removed (ms) 
Derivatives needed to 'whiten' 



52327 - 52538 
2002 Feb 22 - 2002 Sep 21 
78 
52345.0 
15.2238557657(5) 
-4.49522(2) 
2.00(3) 
1.56 
6 



Second phase-coherent solution 



Dates (Modified Julian Day) 
Dates 

Number of TOAs 

Epoch (Modified Julian Day) 

v (Hz) 

v (10- 11 s- 2 ) 



v (10 



-21 



RMS residuals with v removed (ms) 
Derivatives needed to 'whiten' 



52571 - 52776 
2002 Oct 24 - 2003 May 17 
33 
52345 
15.22386798(2) 
-4.5415(1) 
12.33(4) 
1.13 
5 



Third phase-coherent solution 



Dates (Modified Julian Day) 
Dates 

Number of TOAs 
Epoch (Modified Julian Day) 
v (Hz) 



-n 



2 ) 



-21 s -3, 



v (10" 
v (10- 

RMS residuals with v removed (ms) 
Derivatives needed to 'whiten' 



53063 - 54669 
2004 Feb 28 - 2008 Jul 22 
379 
54114.46 
15.21701089718(2) 
-4.48652358(9) 
5.85153(5) 
965 
>12 



a Figures in parentheses are uncertainties in the last digits quoted 
and are the formal la uncertainties reported by TEMPO. 
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Table 3 

Radio - High energy phase offset for 'main' pulse 3 



Pulsar 


Period 


Radio/X-ray 


Radio/Gamma-ray 


Refs. 




(ms) 


onset 


onset 




B1937+21 


1.56 


0.04(1) 


- 


1 


J0218+4232 


2.32 


0.0(1) 


0.50(5) 


2, 3 


B1821-24 


3.05 


0.00(2) 


0.00(5) 


4, 5 


J0613-0200 


3.06 


- 


0.42(5) 


3 


J1614-2230 


3.15 


- 


0.20(5) 


3 


J0751+1807 


3.48 


- 


0.42(5) 


3 


J1744-1134 


4.08 


- 


-0.15(5) 


3 


J0030+0451 


4.87 


0.00(6) c 


0.15(1) 


6, 3 


J2124-3358 


4.93 


— 


-0.15(5) 


3 


J0437-4715 


5.76 


0.003(3) c 


0.45(5) 


7,3 


J0737-3039A 


22.7 


0.04(6) 


— 


8 


Crab 


33.1 


-0.0102(12) d 


-0.001(2) 


9, 10 


B1951+32 


39.5 


— 


0.16(6) 


11 


J2229+6114 


51.6 


— 


0.50(3) 


5 


J0205+6449 


65.7 


0.10(1) 


0.08(2) 


This work, 12 


J1420-6048 


68.2 


-0.35(6) 




13 


Vela 


89.3 


0.12(1) 


0.1339(7), 0.130(1) 


14, 10, 15 


J1028-5819 


91.4 




0.200(3) 


16 


B1706-44 


102 


0.0(2) 


0.211(7) 


17, 10 


J2021+3651 


104 




0.165(10), 0.162(14) 


18, 19 


B1509-58 


151 


0.27(1) 


0.38(3),0.30(5) 


4, 20, 5 


B1055-52 


197 


-0.20(5) c 


-0.25(4) 


21, 22 


B1929+10 


227 


0.06(2) 




23 


B0950+08 


253 


0.25(11) 




24 


B0656+14 


385 


-0.25(5) c 


0.26(8) 


21, 25 


J1119-6127 


408 


0.006(6) c 




26 


B0628-28 


1244 


0.20(5) 




27 



1. ICusumano et al.l J2004T). 2. iKuiper et al.H2004D, 3. lAbdo et alj (12009311 . 4 iRots et al.| 11199811 , 5. IPellizzoni et al l (l2009b|1 , 
6. lAbdo et alJ d2009fl) 7. IZavlin et al.l d2002fl. 8. IChatteriee et al.l [1200711. 9. [Rots et alJ [bo04H. 10. IPellizzoni et al.l (I2009al1. 
11. iRamanamurthy et alj dl995tl. 12. [Abdo et all j2009bl1 13. iRoberts et al.l I2OOIH 14. lHarding et al.l d2002f). 15. lAbdo et all 
ll2009dH. 1 6. lAbdo et al.l fco09cJl. 17 iGotthelf et all ||2002|1, 18. iHalpern et al.l J200SD, 19. lAbdo et all ll2009eD, 20. 
Kuiper et all dl999T>. 21 . iDe Luca et alJ j2005l1, 22 iThompson et al.l 11199911. 2 3. iBecker et~all j200d1. 2 4. IZavlin fc Pavlovl 
25. iHermsen et alj IU997ll , 26. iGonzalez et al.l i2005l) , 27. IBecker et alj ||2005| 1. 



a The phase offset is presented from -0.5 to 0.5 in phase, where a positive value of the phase 
offset indicates that the radio pulse leads the high-energy pulse in phase. Where uncertainties 
are not specified in the original publication, an estimate is given. 

b X-ray pulse represents magnetospheric emission except where otherwise noted. 

c Thermal X-ray pulsations. 

d The phase lag for the Crab pulsar is calculated from the phase of the main peak of the 
radio profile. However, th e Crab pulsar also has a p recursor to the main radio pulse seen only 
at low radio frequencies (Moffctt & Hankins 1996). If the precursor pulse is the true radio 
pulse, the phase offset betwe en radio and the higher energy pulses is ~0.05 in phase (see e.g. 
iRomani fc YadiEaroglulll995l1 . 
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Fig. 7. — Frequency measurements with the post- 
glitch trend subtracted, showing the first observed 
glitch occurring between MJDs 52538 and 52571 
(indicated by vertical lines) with fractional magni- 
tude of Av/v = (3.4 ± 1.1) x 1CT 7 . 1 g uncertain- 
ties in the pre- and post-glitch slopes are shown 
with hatched lines. The change in v as measured 
from the difference in the slope before and after 
the glitch is not statistically significant. 



Fig. 8. — Frequency measurements with the post- 
glitch trend subtracted, showing the frequency in- 
crease occurring between MJDs 52777 and 53062 
(indicated by vertical lines). This could be the re- 
sult of two or more separate glitches occurring dur- 
ing 287-day gap in the data, however, we interpret 
the frequency jump as the likely result of a single 
glitch of fractional magnitude Av/v = (3.8±0.4) x 
10~ 6 . One-CT uncertainties on pre- and post-glitch 
slopes are shown with hatched lines and are in- 
dicative of a change in Av jv = 0.012 ± 0.001. 
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Fig. 9.— Pulse frequency measurements and fitted 
glitch model for the large glitch occurring between 
MJDs 52777 and 53062 for PSR J0205+6449. The 
top panel shows frequency measurements with pre- 
glitch trend removed, with the fitted model for a 
typical glitch epoch over-plotted, while the bot- 
tom panel shows the residuals. The uncertainties 
on each point are dominated by the uncertainty 
in the pre-glitch trend subtracted from the data 
and therefore are of roughly uniform size. For this 
fit to the data, the glitch epoch is MJD 53062, 
with Au/v ~ 3.6 x 10~ 6 , Q ~ 0.66, r d ~ 280 and 

Avjv 0.00453, resulting in X 2 „ = 3-9 for 25 

degrees of freedom. 




Time (MJD) 



Fig. 10. — Pulse frequency measurements and 
fitted glitch model for a typical glitch epoch for 
PSR J0205+6449, considering only the period of 
glitch recovery and ignoring the long-term post- 
glitch change in v. The top panel shows frequency 
measurements with the pre-glitch trend removed, 
with the fitted model over-plotted, while the bot- 
tom panel shows the residuals from the fit. The 
uncertainties on each point are dominated by the 
uncertainty in the pre-glitch trend subtracted from 
the data and therefore are of roughly uniform size. 
Shown here is the best fit model for a glitch oc- 
curring on MJD 52777 with Avjv ~ 5.3 x 10~ 6 , 
r d ~ 383, and Q = 0.67. The reduced X 2 for the 
fit is 0.262 for 12 degrees of freedom. 
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Fig. 11.— Pulse profiles of PSR J0205+6449 for 
RXTE Cycles 6, 7, 9 and 10 for the energy range 
2 - 18keV. Note that RXTE Cycle 6 comprised 
significantly less observing time, resulting in larger 
uncertainties. The interpulse is not clearly visible 
in the Cycle 6 data, but a x 2 test shows that the 
profile is not significantly different from the Cycle 
7, 9, or 10 pulse profiles. 
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Fig. 12.— Pulse profile of PSR J0205+6449 
shown in energy bands 2-10 keV, 10 - 18 keV, and 
18 - 40kcV. Each pulse profile is created by align- 
ing and summing individual profiles from RXTE 
Cycles 7, 9, and 10. Visible structure in the off- 
pulse region of the 18 - 40keV profile is not sta- 
tistically significant as shown by a x 2 test, a Z 2 
test, and a H test. 
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